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ABSTRACT 

The relativistic precession model for quasi periodic oscillations, QPOs, in low mass X-ray 
binaries is reviewed. The behaviour of three simultaneous types of QPOs is well matched in terms 
of the fundamental frequencies for geodesic motion in the strong field of the accreting compact 
object for reasonable star masses and spin frequencies. The model ascribes the higher frequency 
kHz QPOs, the lower frequency kHz QPOs and the horizontal branch oscillations to the Keplerian, 
periastron precession and nodal precession frequencies of matter inhomogeneities orbiting close to 
the inner edge of the accretion disk. The remarkable correlation between the centroid frequency 
of QPOs in both neutron star and black hole candidate low mass X-ray binaries is very well fit 
by the model. QPOs from low mass X-ray binaries might provide an unprecedented laboratory to 
test general relativity in the strong field regime. 
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1. INTRODUCTION 

Old accreting neutron stars, NSs, in low mass X-ray binaries, LMXRBs, display a 
complex variety of quasi-periodic oscillation, QPO, modes in their X-ray flux. The 
low frequency QPOs (~ 1 — 100 Hz) that were discovered and studied from high 
luminosity Z-sources in the eighties are further classified into horizontal, normal and 
flaring branch oscillations (HBOs, NBOs and FBOs, respectively), depending on the 
simultaneous position occupied by a source in the X-ray colour-colour diagram (for 
a review see van der Klis 1995). The kHz QPOs (~ 0.2 to ~ 1.3 kHz) that were 
revealed and investigated with RXTE in a number of NS LMXRBs (see van der Klis 
1998, 1999, 2000 and references therein) involve timescales similar to the dynamical 
timescales close to the NS. A common phenomenon is the presence of a pair of kHz 
QPOs (centroid frequencies of v\ and v<i) which drift in frequency while mantaining 
their frequency difference Av = v<i — v\ xs 250 — 360 Hz roughly constant. Detailed 
studies showed that in four sources Av decreases significantly (by up to ~ 100 Hz) 
as V2 increases; these are Sco X-l (van der Klis et al. 1997), 4U1608-52 (Mendez et 
al. 1998a,b), 4U1735-44 (Ford et al. 1998) and 4U1728-34 (Mendez & van der Klis 
1999). Owing to poor statistics, a similar variation of Av in other sources would 



have remained undetected (Psaltis et al. 1998). 

kHz QPOs show remarkably similar properties across NS LMXRBs of the Z and 
Atoll groups, the luminosity of which differs by a factor of ~ 10 on average. During 
type I bursts from six Atoll sources, a nearly coherent signal at a frequency of 
Vburst ~ 290 — 580 Hz has also been detected (for a review see Strohmayer 2000). In 
a few cases v^ urst is consistent, to within the errors, with the frequency separation 
of the kHz QPO pair Av or twice its value 2Aia Yet there are currently two sources 
(4U1636-53, Mendcz ct al. 1999, and 4U1728-34, Mendez & van dcr Klis 1999) for 
which Vburst is significantly different from Ais and its harmonics. 

The presence of HBOs has been firmly established in both Atoll and Z-sources. 
Their frequency, vhbo (~ 15 to ~ 60 Hz) shows an approximately quadratic de- 
pendence (~ v 2 ) on the higher kHz QPO frequency that is observed simultaneously 
in a number of sources. The frequency changes of the kHz QPOs and HBOs are 
positively correlated with the instantaneous accretion rate. Some evidence has also 
been found for an equivalent of the NBOs and FBOs of Z-sources (Wijnands et al. 
1999; Psaltis, Belloni & van der Klis 1999). 

A remarkable correlation between the centroid frequencies of QPOs (or peaked 
noise components) from LMXRBs has been recently discovered (Psaltis, Belloni & 
van der Klis 1999). This correlation extends over nearly 3 decades in frequency and 
encompasses both NS and black hole candidate, BHC, systems. The frequencies 
of these QPOs, despite their quasi-periodic nature, provide the most accurately 
measured observables of LMXRBs. A primary goal of any QPO model is therefore 
to explain the frequency range and dependence of the different QPO types of these 
sources. The basic features of the relativistic precession model, RPM, are reviewed 
here (Stella & Vietri 1998a, 1999; Stella, Vietri & Morsink 1999). In the RPM the 
QPO signals arise from the fundamental frequencies of the motion of matter in the 
vicinity of the NS. The corresponding orbits are supposed to be slightly eccentric 
and tilted. As in other models, the higher frequency kHz QPOs at v-i are produced 
by the 0- motion (i.e. the Keplerian motion) of inhomogeneities orbiting the inner 
disk boundary, while the lower frequency QPO signal at v\ originates from their 
periastron precession, which is primarily determined by strong-field effects. The 
HBOs are due to the nodal precession in the orbits of the same inhomogeneities, an 
effect which is dominated by frame-dragging around fast-rotating collapsed stars. 
The RPM can be applied to BHCs as well. 

2. PERIASTRON PRECESSION AND KHZ QPOS 

We consider here only infinitesimally eccentric and tilted orbits, under the assump- 
tion that the motion of matter in the innermost disk regions is dictated by the star's 
gravity alone. In the case of a circular geodesic in the equatorial plane (9 = ir/2) 
of a Kerr black hole of mass M and specific angular momentum a, the coordinate 
frequency measured by a static observer at infinity is 

^ = ±M 1 / 2 r- 3 / 2 [2 7 r(l ± aM 1 ' 2 ^' 2 )]- 1 (1) 
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(we use units such that G = c = 1). The upper sign refers to prograde orbits. If we 
slightly perturb a circular orbit in the r and directions, the coordinate frequencies 
of the small amplitude oscillations within the plane (the epicyclic frequency v r ) and 
in the perpendicular direction (the vertical frequency vg) are given by (see Stella & 
Vietri 1999 and references therein) 

u 2 r = 1^1(1 - dMr- 1 ± SaM 1 ^- 3 / 2 -3a 2 r- 2 ) , (2) 

v 2 = v\{\ =f 4aM 1/2 r~ 3/2 + 3a 2 r~ 2 ) . (3) 

In the Schwarzschild limit (a = 0) vg coincides with v^, such that the nodal preces- 
sion frequency u no d = v$ — v$ is identically zero. u r , instead, is always lower than 
the other two frequencies, reaching a maximum for r — 8M and going to zero at 
r ms = QM. This qualitative behaviour of v r is preserved in the Kerr field (a ^ 0). 
Therefore the periastron precession frequency v per = — v r is dominated by a 
"Schwarzschild" term over a wide range of parameters (Stella & Vietri 1999). 

In the RPM the higher and lower frequency kHz QPOs are identified with v 2 = 
v$ and V\ = Vper, respectively. Therefore Ais = v 2 — V\ = v $ — — v T ) = v r . For 
a = 0, Eqs. 1-2 give 

!/ r = ^(l-6M/r) 1 / 2 = ^[l-6(2 7 r^M) 2 / 3 ] 1 / 2 . (4) 

The curves in Fig. 1A show v r vs. v$ for a — and selected values of M, the only free 
parameter in Eq. 4. The measured Av vs. v 2 for eleven NS LMXRBs is also plotted. 
It is apparent that for NS masses in the 2 M Q range, the simple model outlined above 
is in qualitative agreement with the measured values, including the decrease of Ais 
for increasing v 2 seen in Sco X-l, 4U1608-52, 4U1735-44 and 4U1728-34. The model 
above is only an approximation: first, the spacetime around a fast rotating NS is 
different from a Kerr spacetime (due to the star's oblatencss induced by rotation); 
second, the orbits might possess a finite (though small !) eccentricity. Analytical 
formulae to partly correct for these effects were derived by Stella & Vietri (1999); 
Fig. IB shows the fit to the observed Ais vs. v 2 relationship in Sco X-l that was 
obtained through them. The orbital eccentricity, in particular, was varied in order to 
obtain different frequencies, while keeping the periastron distance r p = a(l— e) fixed. 
The best model for a non-rotating NS is shown in Fig. IB. The model reproduces 
fairly accurately the data with a minimum number of free parameters, the NS mass 
(M ~ 1.9 M©) and periastron distance (r p ~ 6.2 M). The latter value is close to 
the marginally stable orbit radius. When the NS spin is allowed a finite value (say 
v s ~ 300 — 600 Hz), fits of very similar quality are obtained, the parameters of 
which differ only slightly from those given above. In essence, the effects induced by 
the NS rotation on v T are small, though non-negligible. 

The behaviour of the curves in Fig. 1 A,B , and therefore the ability of the model 
to match the observations, reflects the properties of the strong field Schwarzschild 
metric, since lower order expansions fail to reproduce the observed frequencies (see 
Stella & Vietri 1999). 
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FIGURE 1. (A) kHz QPO frequency difference Av versus higher QPO frequency 
V2 for eleven LMXRBs. Error bars are not plotted for the sake of clarity. The curves 
give the r- and <fr- frequencies of matter in nearly circular orbit around a non-rotating 
neutron star, of mass 2.2, 2.0 and 1.8 M©. (B) Av versus v-i in Sco X-l. The best 
fit model corresponds to the r- and (^-frequencies of matter orbiting a non-rotating 
1.90 M© neutron star at a periastron distance of 6.25 M (17.5 km). The line marked 
with e gives the orbital eccentricity (xlOOO) as a function of the 4>- frequency. (C) 
kHz QPO frequency difference Av and (double-branched) HBO frequency versus 
higher QPO frequency v 2 in 4U1728-34 (Strohmayer et al. 1996; Ford & van der 
Klis 1998; Mendez & van der Klis 1999). The solid lines give the r-frequency and 
the 2nd and 4th harmonics of the nodal precession frequency v no d as a function of 
the (/(-frequency for infinitesimally eccentric and tilted orbits in the spacetime of a 
1.93 M Q neutron star spinning at 364 Hz (EOS AU; Wiringa et al. 1988). 
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Within the RPM, the maximum value of v r = Ai> depends mainly on the mass of 
the compact object. The NS masses deduced from the simple modelling in Fig. 1A,B 
arc in the ~ 1.8 — 2.0 Mq range, in agreement with the only relatively accurate mass 
measurement from optical spectro-photometry in any of these systems (Cyg X- 
2; M = 1.78 ± 0.23 M Q ; Orosz & Kuulkers 1999). In general within the RPM, 
Av should not be obviously related to the NS spin frequency v s . Therefore, it 
seems natural to identify v s with Vkursti ^.e. the stable frequency seen during type 
I X-ray bursts (for a review see Strohmayer 2000). The distribution of NS spins 
inferred in this way is fairly wide (<~ 290 — 580 Hz) and compares well with that of 
millisecond radio pulsars, MSPs, in agreement with evolutionary scenarios in which 
LMXRBs are the progenitors of MSPs. The 401 Hz spin of SAX J1808.4-3658, 
the only bursting LMXRB displaying coherent pulsations in its persistent emission 
(Wijnands & van der Klis 1998; Chakrabarty & Morgan 1998), is also in the range of 
spin frequencies deduced from Vburst ■ None of the LMXRBs of the (high luminosity) 
Z-class has yet displayed burst oscillations; therefore their spin period is still to be 
measured. Cyg X-2 and GX 17+2, the only type I X-ray bursters in the group, 
might provide this important piece of information. 

3. NODAL PRECESSION AND HBOS 

If the orbits giving rise to the signals at v$ and v per are slightly tilted relative to 
the equatorial plane, nodal precession will take place around the spin axis. In the 
RPM the HBO frequency is related to the nodal precession frequency. From Eqs. 1 
and 3 this can be written in the slow rotation limit (a/M <C 1) 

v nod ~ Anavl ~ 6.2 x lQ- 5 (a/M)mvl Hz ~ 4.4 x 10~ 8 U- > m- x v\v s Hz , (5) 

where M — m M Q . This is the well known Lense-Thirring nodal precession for- 
mula. The latter equality refers to a rotating NS, where aM = 2ttv s I with / = 
10 45 /45 g cm 2 hs moment of inertia. 

If v$ and v s are measured, the only parameter in Eq. 5 that is not identified from 
observations is /45m -1 ; this can vary over a limited range, 0.5 < /45m -1 < 2, for 
virtually any mass and EOS (see the rotating NS models of Friedman, Ipser & Parker 
1986 and Cook, Shapiro & Tcukolsky 1992). The stellar oblateness induced by the 
star's rotation gives rise to correction terms in the nodal precession frequency also 
(see Morsink & Stella 1999 for a post-Newtonian formula). Their relative importance 
increases for high v s and v^. Yet the Lense-Thirring term dominates over a wide 
range of parameters, such that a ~ j/| dependence is expected for v no d- 

An approximately quadratic dependence of vhbo on the higher frequency kHz 
QPOs has been measured in a number of LMXRBs. This dependence was originally 
suggested on the basis of a few power spectra of the Atoll source 4U1728-34 (Stella 
& Vietri 1998a). Ford & van der Klis (1998) analysed a large set of power spectra 
from the same source and determined that the frequency v\ ow of the ~ 10 — 50 Hz 
QPOs scales as ^2' 11±0 ' 06 . Note that the low frequency QPO vs. v<i relation of 
this source appears to be double-branched, with the centroid frequency shifting 
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by a factor of ~ 2 across different observations. This suggests that on occasions 
the 2nd harmonic of uhbo is excited instead of the fundamental. Stella & Vietri 
(1998b) first noticed that the HBO frequency of the Z-source GX 17+2 displays a 
nearly quadratic dependence on v 2 - Psaltis et al. (1999) carried out a systematic 
study of Z-sources and determined that the HBO frequency is consistent with a v\ 
scaling (Cyg X-2 and Sco X-l show evidence for a somewhat flatter dependence). 
In essence these results confirmed one of the basic features of the RPM, namely the 
nearly quadratic dependence of the nodal precession frequency on the (^-frequency. 

If the NS spin frequency is measured, then for any value of the model yields 
a predicted nodal precession frequency which is uncertain only by a factor of a few, 
mainly due to the allowed range of /45m -1 (Stella & Vietri 1998a). Only in the 
Atoll source 4U1728-34 burst oscillations and simultaneous kHz QPOs and HBOs 
have so far been detected unambiguously (Strohmayer et al. 1996; Ford & van der 
Klis 1998; Mendez & van der Klis 1999). Therefore its QPO frequencies can be 
used to test both the vhbo and Av versus v 2 relationships predicted by the RPM, 
when the NS spin derived from burst oscillations is used (vburst — 364 Hz). In order 
to take fully into account of all the effects that contribute determining geodetic 
motion in the vicinity of the NS, we adopted a numerical approach and computed 
the spacetime metric of the star using Stergioulas' (1995) code, an equivalent of 
that of Cook et al.(1992); see also Stergioulas and Friedman (1995). From this, v r 
and v no d were derived as a function of for infinitesimally small tilt angles and 
eccentricities (Morsink & Stella 1999; Stella, Vietri & Morsink 1999). 

Fig. 1C shows the measured values of Av and i^h bo versus v 2 in 4U1728-34. 
Relatively high NS masses (see also Sect. 2) and stiff EOSs such as AU and UU 
(Wiringa et al. 1988) are required in this application of the RPM. The solid lines in 
Fig. 1C are for a 1.93 M© NS with EOS AU and v s = 364 Hz. A good agreement is 
obtained if the HBO frequency, the lower of the two branches seen in 4U1728-34, is 
identified with the 2nd harmonics oiv no d {i.e. 2v no d\ see also Morsink & Stella 1999; 
Stella & Vietri 1999a). Correspondingly the upper HBO branch is well fit by ^v no d- 
The geometry of tilted orbits in the innermost disk regions might be such that a 
stronger signal is produced at the even harmonics of the nodal precession frequency 
(e.g. Psaltis et al. 1999). The frequency range and trend of the epicyclic frequency 
v r in this model are also in reasonable agreement with the Av measurements; a 
more complex model is clearly required in order to fit these data more accurately. 

In summary, the model presented here is capable of reproducing the salient 
features of both the Av versus v 2 and I'm bo versus v 2 relationships, with just two 
free parameters (M and the EOS), the allowed range of which is fairly limited 
(moreover the EOS cannot even be varied continuously !). Concerning Z-sources 
and all other Atoll sources for which burst oscillations have not been detected yet, 
the NS spin can still be regarded as a free parameter. The application of the RPM 
to the HBOs of these sources can therefore be used to constrain the spin of their 
NSs. Conversely, in those Atoll source in which Vburst is measured, but HBOs have 
not been detected yet, the RPM can be used to predict vhbo- These issues are 
briefly addressed in the next Section. 
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4. THE RELATIVISTIC PRECESSION MODEL AND THE PBV CORRELA- 
TION 

Psaltis, Bclloni & van der Klis (1999) recently identified two QPOs and peaked 
noise components the frequency of which follows a tight correlation over nearly three 
decades. This correlation (hereafter PBV correlation) involves both NS and BHC 
LMXRBs spanning different classes and a wide range of luminosities (see the points 
in Fig. 2). In kHz QPO NS systems, these components are the lower frequency kHz 
QPOs, and the low frequency, HBO or HBO-likc QPOs, I'm bo - For BHC systems 
and lower luminosity NS LMXRBs the correlation involves either two QPOs, or a 
QPO and a peaked noise component. In all cases the frequency separation is about 
a decade and an approximate linear relationship (y HBO ~ ^° 95 ) holds. The QPO 
frequencies from the peculiar NS system Cir X-l varies over nearly a decade while 
closely following the PBV correlation and bridging its low and high frequency ends. 
Psaltis, Belloni & van der Klis (1999) noted also that the vs. v\ relations of 
different Atoll and Z-sources line-up with good accuracy. 

The RPM matches precisely the PBV correlation, without resorting to any ad- 
ditional assumption (see Stella, Vietri & Morsink 1999). We assume that in all QPO 
sources, including BHCs, vhbo — 2v no d as in 4U1728-34 (see Sect. 3). Fig. 2A shows 
2v noc i and v$ obtained from Eqs. 1-3 as a function of v per for corotating orbits and 
selected values of M and a/M. The high frequency end of each line is dictated by 
the orbital radius reaching the marginally stable orbit. 

The separation of the lines in Fig. 2A testifies that while v no ^ depends weakly 
on the mass and more strongly on a/M, the opposite is true for By taking the 
weak field (M/r <C 1) and slow rotation (a/M <C 1) limit of Eqs. 1-3 the relevant 
first order dependence is made explicit, 

v+ ~ (2^)- 2 / 5 3- 3 / 5 M- 2 / 5 ^ ~ 33 m- 2 /\ 3 / r 5 Hz , (6) 

v nod ~ (2/3) 6 / 5 7r 1 ' /5 (a/M)M 1 / 5 ^/5 ~ 6 . 7 x iq-2 ( a /M)m 1 / 5 ^/J Hz . (7) 

For the case of rotating NSs we adopt the numerical approach outlined in Sect. 3. 
Results are shown in Fig. 2B for a NS mass of 1.95 M©, EOS AU and v s = 300, 600, 
900 and 1200 Hz (corresponding to a/M = 0.11, 0.22, 0.34 and 0.47, respectively). 
Note that the approximate scalings in Eqs. 6-7 remain valid over a wide range of 
frequencies. Only for the largest values of v per and v s , v no d departs substantially 

6 /5 

from the ~ v per dependence. 

The measured QPO and peaked noise frequencies giving rise to the PBV correla- 
tion are also plotted in Fig. 2B. Higher kHz QPO frequencies from NS systems (^2) 
are included (for the sake of clarity NBOs and FBOs were excluded). The agreement 
over the range of frequencies spanned by each kHz QPO NS system should not be 
surprising: together with the accurate matching of the corresponding v\ — i/ 2 rela- 
tionship in Z-sources, this is indeed part of the evidence on which the RPM model 
was proposed. However the fact that the dependence of v no d on v pev matches the 
observed vhbo ™ v \ correlation to a good accuracy over ~ 3 decades in frequency 
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(down to v\ of a few Hz), encompassing both NS and BHC systems, provides addi- 
tional independent evidence in favor of the RPM. The observed variation of vhbo 
and v\ in individual sources (Cir X-l is the most striking example, see Fig. 2B) 
further supports the scaling predicted by the RPM. The matching of the observed 
v-i vs. v\ relation in terms of v$ vs. v per is also quite accurate. 

For EOS AU and m = 1.95, the v H bo vs - v \ values of most NS LMXRBs are 
best matched for v s in the ~ 600 to 900 Hz range. It is apparent from Fig. 2B 
that v s as low as ~ 300 Hz are required for the Atoll sources with v H bo somewhat 
below the main PBV correlation: these are 4U1728-34 (see Sect. 3) and 4U1608-52 
(from which burst oscillations have not been detected yet). The values above are 
close to the range of v s inferred from Vburst in a number of other Atoll sources 
(van dcr Klis 1999). Z-type LMXRBs appear to require v s in the ~ 600 to 900 Hz 
range, a possibility that is still open since for none of these sources there exists yet 
a v s measurement. Note that the upper HBO branch of 4U1728-34 matches well the 
main PBV correlation. In the interpretation of Sect. 3 the lower branch corresponds 
to 2v noc i and the upper branch to 4v no d- One could further speculate that sources 
following the main PBV correlation, Z-sources in particular, are also in the upper 
HBO branch at ^v no d] in this case their v s might be expected in the ~ 300 — 400 Hz 
range. 

For BHC LMXRBs the scatter around the PBV correlation implies values of a/M 
of ~ 0.1-0.3 (sec Fig. 2B). The points from XTE J1550-564, while inconsistent with 
any single value of a/M , might lie along two distinct branches separated by a factor 
of ~ 2 in vhbOi similar to the case of 4U1728-34. In the RPM the high frequency 
BHC QPOs {e.g. the - 300 Hz QPOs of GRO1655-40) are interpreted terms of 
v\ = v-per- This is at variance with the v\ = v noc i interpretation of Cui et al. (1998), 
which requires high values of a/M (~ 0.95 in GRO1655-40), in contrast with BHC 
accretion-driven spinup scenarios (King & Kolb 1999). From Fig. 2A it is apparent 
that the ~ 300 Hz QPOs from GRO 1655-40 lie close to the high frequency end of 
the a/M = 0.1, m = 7 line. Since the mass of the BHC in GRO 1655-40 determined 
through optical observations is ~ 7 M (Shahbaz et al. 1999), we conclude that, 
according to the RPM, v per ~ 300 Hz close to the marginally stable orbit, where by 
definition v$ — v per . This suggests that any additional QPO signal at V2 = v^ would 
be very close to or even blended with the QPO peak at v\. The detection of two 
closeby or even partially overlapping QPO peaks close to ~ 300 Hz in GRO 1655-40 
would therefore provide further evidence in favor of the RPM interpretation. 

5. DISCUSSION 

5.1. Beat Frequency Models 

In the alternative scenario provided by beat frequency models, BFMs, disk inhomo- 
geneities at the magnetospheric boundary (r m ) and the sonic point radius (r s ) are 
accreted at the beat frequency between the local Keplerian frequency, v^, and the 
NS spin frequency, v Sl giving rise to the HBOs {vhbo = v<p{r m ) — v s ) and the lower 
frequency kHz QPOs {v\ = v<p{r s ) — v s ), respectively (Alpar & Shaham 1985; Lamb 
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FIGURE 2. Twice the nodal precession frequency, 2i/ no d, and (^-frequency, v$, vs. 
periastron precession frequency, v per , for black hole candidates of various masses and 
angular momenta (panel A) and rotating neutron star models (EOS AU, m = 1.95) 
with selected spin frequencies (panel B). The measured QPO (or peaked noise) 
frequencies vi, vi and vhbo giving rise to the PBV correlation are also shown 
in panel B for both BHC and NS LMXRBs and in panel A for BHC LMXRBs 
only; errors bars are not plotted (see Psaltis, Belloni & van der Klis 1999 for a 
complete list of references). We included only those cases in which QPOs at v\ were 
unambiguously detected. NBO and FBO frequencies are not plotted. 
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et al. 1985; Miller et al. 1998). The higher frequency kHz QPOs are attributed to 
the Keplerian motion at the sonic point radius (v 2 — v<p(r s )). The frequency sepa- 
ration Av therefore, yields the NS spin frequency, v s . The narrow distribution of v s 
inferred in this way (~ 250 — 360 Hz) is far from the mass shedding limit of any NS 
model and considerably less extended than that of fast MSPs (~ 600 Hz). Accord- 
ingly the NSs of LMXRBs should be equilibrium rotators with a narrow frequency 
range despite their different average mass transfer rates, magnetic field strcnghts 
and evolutionary histories. Moreover, if LMXRBs are the progenitors of MSPs, then 
BFMs would require a different evolutionary path leading to the formation of radio 
pulsars with v s > 400 Hz. 

The fact that in some Atoll LMXRBs v\> urst ~ Av (or ~ 2Ai/) is readily inter- 
preted, because in BFMs v s = Av. Yet, Av does vary and is significantly different 
from Vburst in several sources (see Sect. 1). This in contrast with the expectations 
of simple BFMs. 

Attempts at fitting the PBV correlation within BFMs by using the range of 
spin frequencies inferred from Av fail to produce a power-law like behaviour over 
a sufficiently large range of frequencies (see van der Klis 1999). This is because 
both v\ and vhbo result from the difference of a variable frequency (y^ at r s and 
r m , respectively) and a fixed frequency (v s ). Moreover, BFMs are not applicable to 
BHCs, since the no hair theorem excludes the possibility that an offset magnetic 
field or radiation beam can be stably anchored to the black hole, as required to 
produce the beating with the disk Keplerian frequency. 

5.2. The Relativistic Precession Model 

The RPM naturally explains the frequency range and dependence of the kHz QPOs 
and HBOs in NS LMXRBs, as well as the PBV correlation, which involves both NS 
and BHC systems. The model has a minimum number of free parameters. Predic- 
tions that can be tested through future analyses and/or observations include: 

(a) The frequency difference Av — v r is expected to decrease also for low values 
of v 2 = v<p (see Eq. 4 and Fig. 1). Moreover if the highest v 2 — v<j> frequencies do 
originate from nearly circular orbits (see Fig. 1A), then Av should quickly decrease 
as v 2 increases further and the orbital radius approaches the marginally stable orbit. 

(b) v 2 = v<p is expected to scale Extending the v 2 vs. v\ correla- 
tion in NS systems toward lower frequencies and detecting the signal at v 2 in BHC 
systems would provide important new tests. 

In the RPM the QPO signals are produced at r/M ~ 100 which, 
for individual sources, must decrease for increasing mass accretion rates (as M 
is positively correlated with e.g. v 2 ). The inferred radii range from close to the 
marginally stable (r/M ~ 6) to r/M ~ 30 over the frequency span covered by 
the PBV correlation. Many NS and BHC LMXRBs display two component X-ray 
spectra consisting of a soft thermal component, usually interpreted in terms of 
emission from an optically thick accretion disk, and a harder, often power-law like 
component, likely due to a hot inner disk region. The QPOs might originate at 
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the transition radius between the optically thick disk and the hot inner region as 
a result of occultation or modulated emission by inhomogeneities (Stella, Vietri &: 
Morsink 1999; see also Di Matteo & Psaltis 1999). 

Mechanisms that can induce a finite (though small) eccentricity and tilt in the 
motion of matter in the innermost disk regions are currently under investigation. 
In the case of NSs, some kind of resonance between the star spin and the motion of 
matter inhomogeneities (see e.g. Vietri & Stella 1998) might be responsible for the 
close commensurability of Vburst an d Ar/ (or 2Af) observed in a few sources. These 
issues will be addressed elsewhere. 

A simple test-particle approximation has been adopted so far within the RPM. 
Much needed hydrodynamical calculations are still in their infancy; among other 
things these are hampered by uncertainties concerning the physics of the innermost 
disk regions. Yet we note that in the hydrodynamical approach explored by Psaltis 
& Norman (2000), the test particle frequencies (the same as in the RPM plus an 
additional frequency at v$ + v r ) are selected by the response of an annulus in the 
disk, when this is subject to a wide-band input noise. 

If confirmed, the RPM will provide an unprecedented opportunity to measure 
GR effects in the strong field regime, such as the periastron precession in the vicin- 
ity of the marginally stable orbit and the radial dependence of the Lense-Thirring 
nodal precession frequency. In principle, accurately measured kHz QPO and HBO 
frequencies would yield crucial information on the compact object such as its mass 
and angular momentum {e.g. by solving Eqs. 1-3 for m, a/M and r). Should suit- 
able, additional observables be found, it might become possible to obtain a self- 
consistency check of the RPM, together with tests of GR in the strong field regime. 
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